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Abstract 

An analytic model for predicting the effect of unsteady local surface injection on the flow separating 
from a streamlined body at angle of attack is proposed. The model uses the premise that separation control 
results from enhanced mixing along the shear layer that develops between the main stream and the fluid 
in the underlying recirculation zone [1]. High-Reynolds-number asymptotic methods are used to connect 
the unsteady surface injection to an instability wave propagating on the separating shear layer and then to 
the large-scale coherent structures that produce the increased mixing. The result is a tool that can guide 
the choice of fluidic-actuator parameters to maximize flow-control effectiveness and may also facilitate 
computer-based numerical experiments. 

1 Introduction 

The drive to reduce NO x and CO 2 emissions from current and future commercial aircraft has led to an increased 
emphasis on improving the efficiency of turbofan engines. The performance of modern inlet/fan/compressor systems 
has for the most part reached the limits achievable by geometric shaping, i.e. by simply bending metal. To realize 
further gains, losses that occur locally within the turbomachinery blading and/or transiently during certain operating 
ranges must be addressed. 

One common source of loss is low-momentum regions that develop along the suction sides of the aerodynamic 
surfaces which makeup the various engine components. In the compressor, the low-momentum regions give rise to a 
stall boundary which limits the amount of work obtained per stage and so dictates the number of stages required to 
achieve a desired pressure rise. The ability to manage this loss through adaptive flow control has now become possible 
due to recent advances in micro electro-mechanical systems. 

Numerous experiments on isolated airfoils in cross flow and a smaller number of investigations of flow through 
blade cascades have demonstrated the potential of local surface injection to delay or prevent flow separation. Com- 
plementary studies using unsteady injection have shown that the same measure of flow control can be achieved with 
significantly less injected mass. 

Limiting injected mass and thereby thermodynamic cycle penalties incurred by bleed extraction is just one of 
several objectives which must be considered when designing intelligent flow-control approaches for aeropropulsion 
systems. Additional considerations include minimizing the net injected momentum to avoid parasitic mixing losses 
and optimizing the acoustic characteristics of the fluidic actuator to reduce power consumption and improve frequency 
response. 
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Designing a robust effective flow-control system that maintains optimum efficiency throughout the engine oper- 
ating envelope requires tailoring the actuation device for the specific environment in which it is placed. Unfortunately, 
parametric studies of pulsed or time-harmonic injection systems are difficult to perform experimentally (even in single 
blade configurations) because the useful operating range of any particular actuator is dictated by the geometries of its 
attendant chamber and orifice. Time-accurate numerical simulations are only marginally less difficult since altering 
the actuator design generally requires re-griding the problem. Similarly, the low-order models for the time-averaged 
effect of unsteady surface injection required by existing multi-stage steady-flow design codes are lacking. 

The present investigation addresses these issues by providing an analytic model that predicts the effectiveness of 
separation control as a function of tunable actuator parameters. Asymptotic methods are used to formalize a mecha- 
nism for separation control proposed in reference [ 1 ] wherein the basis is mixing enhancement along the separating 
shear layer. The performance of a given flow-control device is measured in terms of the ratio of mass entrained by 
the separating shear layer to mass available in the recirculation zone. The analysis is done in the framework of the 
high-Reynolds-number limit which is a reasonable approximation for the flows of interest here and, in any event, often 
correctly captures the physical balances in play at moderate Reynolds number. 

The problem is formulated in §2 where the asymptotic solution for the time-mean flow that would exist in the 
absence of control is given near its separation point. The unsteady flow field is analyzed in §3 where the coupling 
between the component forced by the surface injection and the instability wave excited on the separating shear layer is 
considered. In §4, a model for the enhanced mean-flow spreading rate along the shear layer is given and a flow-control 
metric is proposed. Some preliminary results obtained from the model are presented and their implication discussed 
in 85. 


2 Formulation 

At frequencies typical of turbomachinery applications, unsteady surface injection is thought to effect separation 
control through mixing enhancement. Surface injection introduces an unsteady perturbation to the flow that excites 
a spatially growing disturbance on the separating shear layer by the Kelvin-Helmholtz instability mechanism. The 
instability wave grows as it propagates downstream until nonlinear effects cause it to roll up into a discrete vortex. 
Further downstream, the individual vortices, which are identified with large-scale coherent structures, begin to merge 
and fine-grained turbulence becomes the dominant vehicle for shear-layer mixing. This evolutionary process is ac- 
companied by enhanced entrainment of low-momentum near- wall fluid which tends to drive the separating shear layer 
toward the surface. If the rate of entrainment is sufficient, the streaming flow re-attaches. 

The above description is cast into a workable tool for flow-control analysis by combining existing high-Reynolds- 
number analytic models for the following distinct fluid-mechanical phenomena that dominate various stages of the 
flow development: 1 ) the uncontrolled base flow, i.e. the time-averaged separated flow to which control is applied; 2 ) 
the linear unsteady flow field which includes the perturbation associated with the imposed surface injection and the 
shear-layer instability which it excites; 3) the large-scale coherent mode into which the instability evolves; and 4) the 
enhanced mixing produced by the coherent mode. 

2.1 Time-Averaged Base Flow 

Consider the two-dimensional streaming flow of a compressible fluid past an airfoil of length L. The pertinent 
situation arises at angles of attack at which (absent intervention) the flow breaks away from the suction side of the 
airfoil an order L distance downstream of the leading edge. Let U s , R s , P s and n s denote the characteristic speed, 
density, static pressure and dynamic viscosity along the separating streamline, respectively. A subsonic adiabatic 
flow is assumed and, since the analysis to follow is limited to a small neighborhood of the separating streamline, the 
pressure-density relation P = R 1 is approximated by the Chaplygin-Karman-Tsien equation of state [2, p. 285] 

p ~ l = 1 ^ir (1) 

where pressure P and density R are non-dimensionalized by P s and R s , respectively, and 7 is the ratio of specific 
heats. Equation (1), which is just a tangent approximation to the exact pressure-density relation, greatly simplifies the 
analysis by allowing use of the hodograph method. 

The limiting form of the flow field on the airfoil scale as the Reynolds number Re = R S U S L/ fi s tends to infinity 
can then be found in the class of potential flows described by a generalized version of Kirchhoff’s free- streamline 
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Figure 1 . Separated flow configuration. 


theory [3]. The solutions are parameterized by the point of intersection between the airfoil surface and the free- 
streamline that separates the oncoming flow from the stagnant fluid in the recirculation zone, i.e. regions 1 and 2 of 
figure 1, respectively. References [4] and [5] show how the exact location of this separation point depends on the 
detailed nature of the viscous flow along the airfoil surface. 

The base-flow separation point is left unspecified at present, but it is assumed that both the site of fluid-injection 
and the onset of shear-layer roll-up lie within a short distance (relative to L) of that point. Accordingly, it suffices to 
give the base-flow solution only in the vicinity of its separation point which is now taken as the origin of the Cartesian 
coordinate system {X, Y} = X with X and Y tangent and normal, respectively, to the surface (cf. figure 1) and both 
non-dimensionalized by L. 

The Kirchhoff theory, when generalized to compressible fluids satisfying (1), predicts the shape of the free stream- 
line separating regions 1 and 2 to be 

Y = H S {X) = aX 3/2 + bX 2 + 0(X 5 / 2 ) (2) 

for 0 < X <C 1, where a is a constant that depends on the viscous flow along the surface and b is half the non- 
dimensionalized airfoil curvature at X = 0. References [4] and [5] show that for high-Reynolds-number separations 
from a smooth surface 0 < a ^ 0(1). The separations of interest here result from the action of an adverse pressure 
gradient so b < 0 and, to required order of accuracy, the airfoil shape is 

Y = H b (X) = bX 2 + 0(X 3 ) (3) 

for \X\ <c 1. 

In region 1, the base-flow velocity (non-dimensionalized by U s ) is given as 

d#i 

U 1 ={U 1 ,V 1 } = i+^r, (4) 

with velocity potential 

s ' = -7 t^ 5 r *(' z3,2 ) + °(i z ' 2 » < 5 > 

for \X\ <C 1, where i is the unit vector in the X direction, i = 1, Z = X + i^/l — Y, M s = U s /y/jP s /R s 

is the characteristic Mach number along the free streamline and a numeric subscript on a dependent variable indicates 
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the flow region to which it applies. Corresponding expressions for the region 1 density and pressure can be found from 
(1) and the appropriate Bernoulli equation, eg. 

-Pi = 1 + 7-V? (1 — U\) + 0(\Z\). 

In region 2, the velocity vanishes and the density and pressure take on their constant free- streamline values. 


3 Linear Unsteady Flow Field 

The unsteady fluid motion generated by the surface injection (but external to the actuation device) introduces 
an additional length scale U s /ujj into the problem where ujj is a characteristic angular injection frequency. In flow- 
control theory, it is usual to work in terms of a reduced frequency (or Strouhal number), eg. F + = ujjL/2irU s . 
Values of reduced frequency reported in the literature vary widely however, for flow-control approaches where mixing 
enhancement is the active mechanism, F + is approximately order one. (At higher reduced frequencies, effective flow 
control is thought to result from virtual aeroshaping.) In the following analysis 


5(Re) < e = 


U s j u) j 

L 


1 


27rF+ 


< 1 


( 6 ) 


is assumed where the lower limit is introduced so that viscous effects are for the most part passive. The precise form of 
T depends on the detailed nature of the viscous flow along the airfoil surface, eg. T = i?e -1 / 4 in the case of laminar 
separation. 

Minimization of parasitic losses is a primary concern in the design of efficient flow-control systems. With this 
in mind, the injection strength (measured in terms of net momentum flux, say) is assumed sufficiently small for 
the unsteady dynamics to remain predominantly linear at least within an order U s /ujj neighborhood of the base-flow 
separation point. The unsteady-flow field is then characterized by time scale t non-dimensionalized by Wj and includes 
a spatial variation on the scale x = X/e which is short compared to the airfoil chord. 

It follows from these considerations that the velocity field (external to the actuation device and non-dimensionalized 
by U s ) can be expressed as 


u i — SilUl + V 


fa + SuMg 




(7) 


for i = 1, 2, with potentials fa governed by 






( 8 ) 


where Sij is the Kronecker delta function, V = d/dx and similar expansions apply for the density and pressure. 
Equation (8) must be solved subject to slip or injection conditions on the airfoil surface ey = as well as appropriate 
far-field conditions as y — » oo. 

Solutions in regions 1 and 2 are, of course, coupled by kinematic and dynamic conditions at the free streamline. 
Upon eliminating the unsteady perturbation to H s and expressing the result in terms of ( fa these conditions become, to 
required order of accuracy, 


d \d ± ( d d \ 

^[di^~\di + di) 


d 

at 


d d 
dt + dx 


* i ~ I* 


= Lk (4> I,<h; eH'J, M 2 S ) , 


(9) 

( 10 ) 


at ey = H s for x > 0, where L K and L D are linear operators on fa (see appendix) and a prime denotes differentiation 
with respect to X. 

The linear nature of the unsteady problem is now exploited by expressing the solution as the sum of two compo- 
nents: one directly forced by, but external to, the actuator and the other excited along the separating shear layer by 
the Kelvin-Helmholtz mechanism. Both components are, of course, ultimately related to the surface injection, but the 
(artificial) splitting is used here to simplify the subsequent analysis by isolating the distinct spatial scales on which the 
unsteady flow evolves. 
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3.1 Forced Solution 

Reference [6] shows that small-amplitude unsteady injection along a flat plate over which a high-Reynolds number 
flow passes gives rise to both acoustic and hydrodynamic modes of motion. The latter is confined to the thin viscous 
layer adjacent to the surface and takes the form of Tollmien-Schlichting waves. The T-S waves will grow as they 
propagate downstream if the injection frequency lies within the unstable band dictated by the local boundary-layer 
profile. The assumption employed here and embodied in (6) and (7) is that only the acoustic motion persists at order 
Us/uj distances from the injection site, i.e. that c Oj is less than the so-called lower branch frequency. 

Attention is focussed on injections through two-dimensional slots of width dj placed a distance £j upstream of 
the base-flow separation point (cf. figure 1) where typical values encountered in synthetic and pulsed jet flow-control 
experiments suggest dj <C £j with £j order U s /ujj. For simplicity, the injection is assumed to be normal to the airfoil 
surface with a monochromatic velocity of period 2tt / uj. More complicated injections can always be obtained by 
superposition. 

The leading-order outward propagating solution to (8) in region 1 is then 


<t>{ ~ 


1 1 djLOj 

8 ^/T^mI~u7 



(t + M? 


X — Xj \ 


Hi, 1 ’ 



1 ~Mf) 


+ c.c. 


( 11 ) 


where superscript / indicates the forced component, Vj (t) is the area-averaged slot exit velocity non-dimensionalized 
by U s , H£ 1} is the Hankel function of the first kind and zero-th order, z = Z/e and Xj = —£jWj /U s . 

The right-hand sides of (9) and (10) give rise to an order e 1 / 2 correction to (11) which, owing to the behavior of 
H”, becomes singular as \x\ — > 0. It is the removal of this singularity by application of an unsteady Kutta condition 
that provides the coupling between the forced and excited components of motion. 


3.2 Excited Solution 

The component of unsteady fluid motion excited along the separating shear layer is more strongly dependent on 
the flow in the recirculation zone than the forced component. It is only when the transverse extent of region 2 becomes 
order U s /ujj that the mechanism responsible for the excited motion takes on the classic Kelvin-Helmholtz instability 
form. Equations (2) and (3) show that this stage is reached when X becomes order (e/a ) 2 / 3 at which point the excited 
solution varies (locally) on the x scale. 

As indicated above, the coupling between the forced and excited solutions takes place near the base-flow sep- 
aration point as | a? | — ► 0. It is therefore necessary to determine the excited solution over the stream wise scale 
x\ = (a/e) 2 / 3 X where the curvature of the base-flow free streamline becomes important. If the dynamics remain 
linear on the x\ scale, the uniformly valid solution to (8) in region 1 has the form 


^1 rSJ A(xi) exp 


i(a 2 e) -1 / 3 


nx i 

/ a(£)d£ ~ Ki(y~ e-'Hg) - it 

Jo 


c.c. 


( 12 ) 


for x\ > 0, where superscript e indicates the excited component, the slowly varying amplitude function A is de- 
termined by a solvability condition at next order (cf. figure 2) and c.c. denotes the complex conjugate. The slowly 
varying wavenumber a(x \ ) is determined by the well-known characteristic equation describing the Kelvin-Helmholtz 
instability of a compressible vortex sheet near a solid boundary [7] which in the present notation reads 


(1 — a) 2 K 2 tanh 



+ /£ 1 — 0 


(13) 


where Ki = with branch cuts fixed by boundedness and propagation considerations in the 

respective regions. 

The root to (13) of relevance here is the one that yields a spatially growing downstream propagating instability 
wave, i.e. the root with Rea > 0 and Ima < 0 (cf. figure 3). The expansion of the corresponding eigenfunction 
(to which (12) represents the leading-order term) breaks down as x\ becomes small and it is necessary to analyze the 
flow in an inner region where \X | is order (ae 2 ) 2 / 3 and the dynamics are quasi steady. The leading-order inner-region 
solution turns out to have a singularity at the origin of the same form as that exhibited by the forced component. 
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e~ l (H s -H b ) 


Figure 2. Slowly varying instability-wave amplitude as function of scaled region 2 thickness. □, M s = 0.2; A, M s 
M s = 0.8. 



Figure 3. Unstable root to (13) at M s = 0.6 as function of scaled region 2 thickness. Open symbols, Rea; filled symbols 
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The analysis of reference [8] is used here to conclude that viscous (and interactive) effects act locally to smooth out 
the combined singularities of the complete flow (i.e. steady base flow plus forced and excited unsteady components) 
such that when viewed on the comparatively large scale of the inner region the singularities of the forced and excited 
solutions exactly cancel. The result is an unsteady Kutta condition that relates the instability wave amplitude A to the 
surface injection as follows 


where 


3 1 dj L j ( M 2 a L\ 4/3 n 

128 m^/T^mI T Tyi-MiTJ b 


S(w) = w 1 / 3 e -i(M s a.- 7 r/ 3 ) 

and Lj is the stroke length based on the 
from the slot during the discharge stroke. 


1 + i 


. 1 + M; 


M 9 


+ (3 - 


area-averaged injection velocity, i.e. 



i^-i ) 

l-M 2 ) 


(14) 


2 M 2 + M s 4 ) wH^u;) 

the length of the slug of fluid pushed 


4 Mixing Enhancement 

The exponential growth of the excited component of unsteady motion causes it to dominate the forced component 
downstream of separation. This rapid growth continues until nonlinear effects come into play to produce a spatial 
equilibration as the instability wave rolls up into a discrete vortex [9] . It is generally believed that such linear instability 
waves are the progenitors of the large-scale coherent structures which dominate the spatial development of transitional 
and turbulent shear flows. 

Interest here is in the effect injection-induced large-scale coherent structures have on momentum transfer across 
the separating shear layer. A time-averaged measure of the fluid entrained by the shear layer is provided by the mean 
shear-layer momentum thickness S(s) where 8 is the downstream distance along the centerline and both s and S are 
non-dimensionalized by L. The conservation equations for the mean (i.e. Reynolds averaged) motion then show that 
the growth in S caused by viscous diffusion is augmented by Reynolds stresses due to both fine-grained turbulence and 
large-scale coherent structures. 


4.1 Mean-Flow Spreading Rate 

Reference [10] presents a high-Reynolds-number, energy-transfer-based evolutionary model for the mean-flow 
spreading rate d 5 / ds that includes effects due to viscous dissipation and the production of and interaction between 
coherent modes and fine-grained turbulence. The transverse integrals of the Reynolds-stress terms which appear in 
the model are evaluated by exploiting the ability of local linear inviscid stability theory to accurately predict the fast- 
scale propagation and cross-stream distribution of the coherent structures even in the presence of turbulence. The 
corresponding amplitudes are however determined by nonlinear considerations on the s scale. 

The most rapid rise in spreading rate (and hence entrainment) occurs downstream of the linearly unstable re- 
gion where viscous diffusion gives rise to a y / s growth in S, but upstream of the vortex-merging stage where fine- 
grained- turbulence production leads to a linear growth in S. For a monochromatic coherent-mode in the intermediate 
streamwise region, the evolutionary spreading-rate model simplifies to 

i ^ = \A\ 2 I rss I = \A\ 2 I rs , (15) 

ds ds 

where A(s) is a nonlinear amplitude function based on appropriately normalized eigenfunctions from linear stability 
theory, / is the mean-flow kinetic-energy flux integral, I is the time-mean of the kinetic-energy flux integral for the 
coherent mode and \A\ 2 1 rs gives the energy exchange rate between the mean flow and coherent mode. The flux 
integrals / and / are constants while I rs is a function of S [10]. The latter is approximated here as 

I rs « CS - DS 2 

where, in general, C and D depend on s. Equations (15) are integrated to yield 

Soo-S _ - S\A\ 2 

<5oo-<5 0 ^ool^ool 2 - ^o|^4o| 2 ( 
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Figure 4. Normalized flow-control metric as function of scaled downstream distance with £ = 0.006, # = 0.25 and k r X 0 = 9. 


and 


S\A\ 2 


<5q|A)| 2 <5oo|^oo| 2 

(5 0 |^4 0 | 2 + OoolAool 2 - <5o|A)| 2 ) exp ^—2 f s s o Mo) 


(17) 


where k r = S 00 \A 00 \ 2 D/2I is an 8 dependent local linear growth rate, £oo|Ax)| 2 = ^o|^o| 2 + (^oo — Sq)I/I and 
Soq = C/D are equilibration values that depend on both frequency and Reynolds number and subscript 0 indicates 
evaluation at a stream wise location so near the end of the linearly unstable region. 

4.2 Flow-Control Metric 

The degree to which the injection induced mixing enhancement effects separation control is quantified using a 
premise presented in reference [1]. Accordingly, the oncoming flow re-attaches to the airfoil if and when the mass 
of fluid entrained into the shear layer exceeds that available for entrainment, i.e. the fluid in the recirculation zone. 
The amount of mass entrained M e over an arbitrary downstream distance is proportional to the integral of S over the 
corresponding distance along the shear-layer center line while the amount of mass available M a is proportional to the 
corresponding volume of region 2. 

To significantly impact airfoil performance, the re-attachment must occur a short distance (relative to chord) 
downstream of the separation point. Consequently, the small-X approximation is made and s = X + 0(X 4 ). If it is 
further supposed that k r is effectively constant over the region of interest then it follows from (2), (3), (16) and (17) 
that the flow-control measure of merit is 

M e ^ bSpo (1 - g)2 k r (X - Xp) + (!-#) In {e + (1 - g) exp[-2 k r (X - X 0 )}} 

M a ^ dakr (! _ £ )(j5/2 _ x^ /2 ) 1 

where $ = ^o|^o| 2 /^oo|^oo| 2 and $ = ^o/^oo are disturbance-energy and mean shear-layer thickness parameters, 
respectively. 
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Figure 5. Normalized disturbance energy as function of scaled downstream distance with e = 0.006. 


5 Results and Conclusions 

The flow-control metric determined by (18) is plotted in figure 4 in normalized form for values of 6, d and k r X o 
estimated from the experimental data considered in reference [9] . The figure reveals the competition between shear- 
layer entrainment and recirculation-zone expansion as the flow evolves downstream of separation. At small abscissa, 
the shear-flow spreading rate exceeds the rate of recirculation-zone expansion and the flow-control metric increases. 
This trend continues until the instability wave roll-up and spatial equilibration are complete (cf. figure 5). The en- 
hanced mixing generated by the coherent mode is then maximal and further recirculation-zone growth results in a 
decline in the flow-control metric. Accounting for vortex merging and its attendant fine-grained turbulence production 
which are dominant at the larger downstream distances merely alters the rate of decline. Overall the curve in figure 4 
is in good agreement with the experimental data presented in reference [1, fig. 4]. 

Under the premise employed here, successful separation control results when the ratio M e /M a exceeds unity. 
The location of the peak in figure 4 is decisive. In this regard, it should be noted that both abscissa and ordinate 
depend implicitly on the fluid injection parameters. In the simplest approximation, k r is given by the linear weakly 
non-parallel analysis of §3.2, so 


k r (X - X 0 ) » -Ima(oo) y^{X - X 0 ), 

Us 

and the abscissa is proportional to a reduced frequency / + = F + (X — Xo) based on the stream wise length of the 
recirculation zone. The peak in M e /M a is then suggestive of an optimal injection frequency, however the dependence 
on Xo means that the optimum is not global. The normalization (M e /M a ) o used for the ordinate is proportional to 
the ratio of shear-layer to recirculation-zone thickness at X = Xo and thus varies as l/Xo. Shifting the onset of 
nonlinearity upstream (i.e. reducing Xo) increases the likelihood of flow re-attachment. 

The quantity Xo is determined by the properties of the instability wave from which the coherent mode originates. 
It can be reduced by increasing either the initial linear instability wave amplitude or frequency. Relating Xo to 
the injection parameters is, however, complicated by the coupling condition (14) which shows that the initial linear 
amplitude A(0) involves its own frequency dependence. Figure 6 is a plot of |A(0)| as a function of frequency for 
various M s . It follows that Xq decreases with slot width dj and stroke length Lj, increases with distance between 
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Figure 6. Initial instability-wave amplitude as function of injection frequency. □, M s = 0.2; A, M s = 0.6; ©, M s = 0.8. 


injection site and separation point tj and can either increase or decrease with non-dimensional injection frequency 

£j(jjj/U s . 

Clearly, even in the context of the foregoing simplifications, the flow-control metric (18) involves an intricate 
dependence on the fluid injection parameters. This may explain the wide range of reported frequencies at which un- 
steady surface injection successfully delays or prevents flow separation. Understanding this parametric dependence is 
necessary in order to guide the optimization of actuator design for maximal flow-control effectiveness while accurately 
capturing it in a deterministic-stress model is a prerequisite for performing system-level numerical predictions using 
steady-flow design codes. 
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Appendix 

The linear operators on the right-hand sides of (9) and (10) are given as 
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to required order of accuracy. 
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